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We investigate epidemic models with spatial structure based on the cellular automata method. 
The construction of the cellular automata is from the study by Weimar and Boon about the reaction- 
diffusion equations [Phys. Rev. E 49, 1749 (1994)]. Our results show that the spatial epidemic 
models exhibits spontaneous formation of irregular spiral waves at large scales within the domain 
of chaos. Moreover the irregular spiral waves grow stably. The system also shows spatial period- 2 
structure at one dimension outside the domain of chaos. It is interesting that the spatial period-2 
structure will break and transform to spatial synchronous configuration in the domain of chaos. 
Our results confirm that populations embed and disperse more stably in space than they do in 
non-spatial counterparts. 

PACS numbers: 05.50.+q, 87.23.Cc, 87.18.Hf, 89.75.Fb 



I. INTRODUCTION 

Several theoretical models have shown that population 
invasion and dispersion is more stable in space than that 
in non-spatial counterparts. More stable mean that a 
previously unstable equilibrium point becomes stable un- 
der a greater variety of conditions, or that an equilibrium 
is approached faster [1]. In oscillatory systems where the 
equilibrium is a limit cycle or more generally an unstable 
focus, diffusion or dispersal will create wave-like patterns. 
The Lattice Lotka-Voltera (LLV) model was studied ex- 
tensively UMM' 

In non-linear ecology system, the two 
most commonly seen patterns are spiral waves and turbu- 
lence. Spiral waves play an important role in ecological 
systems. For example, spatially induced speciation pre- 
vents extinction for the predator-prey models [H, @] • A 
classical epidemic model is ordinary differential equations 
(DDEs), or called mean- field (MF) approximation [7]. 
The DDEs methods (or MF) are based on the assumption 
that the population is well-mixed, with the subpopula- 
tions involving susceptible, infected, removed, etc., inter- 
acting in proportion to their sizes. Nonspatial theory typ- 
ically predicts selection for maximal number of secondary 
infectors. Among these epidemic features, the existence 
of threshold values is crucial for the spread of an infec- 
tion [8|,^]. A second classical approach describes spatially 
extended subpopulations, such as a coupled map-lattice 
models [To^, reaction-diffusion equations, deterministic 
cellular automata and integrodifference equations model. 
In the literature [11] the author argues the numerical 
simulations in a predator-prey system, and shows that 
there are either irregular spatiotemporal oscillations be- 
hind the invasion or regular spatiotemopral oscillations 
with the form of a periodic travelling 'wake' depending 



on parameter values. 

In our paper, the geographic spread of an epidemic 
can be analyzed as a reaction-diffusion system, in which 
both subpolulation exhibit local random movement, and 
the al gori thm of cellular automata are based on these 
paper jlll . [l2| . More recently, studies have shown large- 
scale spatiotemporal patterns in measles [l3| and dengue 
fever [lj|. These studies have shed new light onto key 
research issues in spatial epidemic dynamics, but the de- 
tailed theoretical studies are difficult. The study of pop- 
ulation dynamics takes into account the species distri- 
bution in space, interactions between individual species 
that are located in the same neighborhood, and mobility 
of the various species jTH, [H, These studies pre- 

dict the formation of spatial complex structure, phase 
transitions, multist ability, oscillatory regions, etc. In 
the Ref. [17], the author studies the susceptible-infected- 
resistant-susceptible (SIRS) models with spatial struc- 
ture using cellular automata rules, showing the forma- 
tion process of the spatial patterns (turbulent waves and 
stable spiral waves) in the two-dimensional space and ex- 
istence of stable spiral waves in the SIRS model. 

The principal objectives of the present work is that 
the susceptible-exposed-infected-resistant (SEIR) model 
with spatial structure is investigated by using cellular 
automata algorithm. The SEIR model and its classical 
DDEs version are presented in [18',T9']. In fact, many dis- 
eases are seasonal, and therefore an important question 
for further studies is how seasonality can influence spa- 
tial epidemic spread and evolution. Hence, we consider 
the seasonal parameter, f]{t) = /^o(l + ^ sin(27rt), where e 
is the fluctuating amplitude of contact rate. Commonly, 
we describe the susceptibility, exposure, infection and re- 
cover process in terms of four nonlinear DDEs. We use 
S for susceptibles, E for the exposed, / for infectors, and 
R for the recovered. The dynamical equations for SEIR 
model are 
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Here /i is the death rate per capacity, 1/S and I/7 are the 
mean latent and infectious periods of the disease. /3{t) is 
the rate of disease transmission between individuals. The 
population can be normalized to S-\-E-\-I-\-R = 1, so all 
dependent variables represent fractions of the population. 
The original studies show that the system of ([T]) exist 
three phase transitions which are the stable behavior, 
the limit cycle and the chaotic behavior in the mean- 
field limit [18|] with respect to the fiuctuating amplitude, 

€. 



II. NEIGHOURHOOD-DEPENDENT MODEL 

Generally, studies on the spacial epidemic models show 
that there exists spatio-temporal travelling waves [20| 
(e.g. dengue haemorrhagic fever (DHF) [13, Hll] and 
measles [13]). However, few systems are well enough doc- 
umented to detect repeated waves and to explain their 
interaction with spatio-temporal variations in population 
structure and demography. The actual epidemic spread 
is spatio-temporal and local individual interact. Here we 
study the individual moving of the susceptible, exposed, 
infector and recover, and their diffusion from one lattice 
site to another. Then the equations ^ read 



dS{r,t) 
dt 

dE{r,t) 
di 

dt 

dR{r,t) 
dt 



= /i - f3{t)IS -fiS^ Di\/^S{r, t). 
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: P{t)IS - (/i + S)E + D2V^E(r, t), (2b) 



= ^^-(7 + /i)/ + I^3V2/(r,t), 
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We study the system (|2]) using cellular automata 
method, which is suitable for modeling many reaction- 
diffusion systems in a quantitatively correct way based 



on the Ref. [12], and demonstrate recurrent epidemic 
spiral waves or traveling waves in an exhaustive spatio- 
temporal through the numerical simulation. Simply, we 
use c(r, t) to denote the vector of individual density in 
position r and at time t, and L{c{r^t)) to describe the 
local kinetics; D is the diffusion coefficient matrix. Then 
the system (|2]) can be written as 



dc{r,t) 
dt 



L{c{r,t))^DV'^c{r,t). 



(3) 



In the following simulation, we may discard the equa- 
tion (|2dp . since we concern the susceptible, exposed and 
infected. 

We define this model as follows. Space is made up of a 
square lattice of J x J. In each step the individuals ran- 
domly move in its neighborhood. The state of the cellular 
automata is given by a regular array of density vector c 
residing on a two-dimensional lattice. We consider cel- 
lular automata with hi states (denoted by the integers 
0, 1,2,3, ••• ,6i). Here species state and hi are zero 
population level and maximum population level respec- 
tively. The first step of each time iteration corresponds 
to local dynamics, and the state at each spatial lattice 
changes independently of the states at other vicinity lat- 
tices. The second part of each time step corresponds to 
unbiased spatial movement. The central operation of the 
cellular automaton consists of calculating the sum 



Ci{r,t) = ^ Ci(r + r'). 



(4) 



r'eNi 



Where the summation takes up all of the nearest neigh- 
bors of the cell r. The neighborhoods can be different 
for each species i. We use the Moore neighborhood for 
all i in the two-dimensional space, i.e. 



Ns,uare = {(0, 0), (1, 0), (0, 1), (-1, 0), (0, -1), 

(1,1), (-1,1), (1,-1), (-1,-1)}. (5) 



We normalize values of Ci(r,t), and the Ci{r^t) = 
Ci(r, t)/{hiNi) is local average density of the Ci(r, t). The 
Ci(r, t) is always between zero and one. 

From the Ref. [l^], the two-dimensional discretization 
version of Eq. (|3|) takes the form 



Ci(r, t + 1) = Q(r, t) + ML{ci{r, t)) + AV^Q(r, t), i = 1, 2, 3, 4. 

I 



(6) 



Where the 



and Di defines the space scale. 



Di=Dii^, ^ = 1,2,3,4. 



(7) 



3 



Furthermore we have 

Ci{r,t+1) = L*{ci{r,t)), 



1,2,3,4. 



(8) 



Where L*{ci{r,t)) = Ci{r,t) + AtL{ci{r,t)). 

As c(r, t) is the average output of the CA for the sys- 
tem and therefore it is given by 



Cjir,t + 1) = 



*( Ci{r,t) \ 
\ hiNi )\ 



(9) 



for species j. About c(r, t) as output of the CA described 
in detail can be found in the paper [l2[ . 



III. NUMERICAL RESULTS 

We have performed extensive numerical simulations of 
the described model, and the qualitatively results are 
shown here. In cellular automata simulation, periodic 
boundary conditions are used and At = 0.005. The space 
scale Di = 0.2, D2 = 0.05, Ds = 0.02 and grid size used 
in the evolutional simulations is 100 x 100 cells. Ev- 
ery species has 100 states in the system (|2]), and more 
states enable more accurately for discrete representation 
of the continuum models, while it is complex for analysis, 
and this is described in detail elsewhere [22]. We have 
tested that the larger grid size does not change qual- 
itative result for the evolutional dynamics. The con- 
tact rate fluctuates with the seasons can be approxi- 
mated in several ways. Simply we choose sinusoidal force, 
P{t) = Po{l + £:sin(27rt), where < 6 < 1, and another 
more realistic option is term-time force, which sets trans- 
mission rates high during school terms and low in other 
place [23I . The spatial patterns evolve from random ini- 
tial conditions. The maximum density of susceptible lev- 
els, exposed levels and infected levels are set to 50, 5 and 
1, respectively, in the two-dimensional space. Other ini- 
tial conditions have been explored as well, and no change 
has been observed in the behavior. In Fig. [1] and Fig. [3] 
four different snapshots during the temporal evolution 
of the system are presented in two-dimensional space. 
These figures and the following figures are species den- 
sity levels as a function of space and time on a gray scale, 
with white corresponding to the lowest-density state and 
black corresponding to the highest density state. 

Fig. [U and Fig. [3] have depicted spatial patterns in the 
two-dimensional space under different e values (the fluc- 
tuating amplitude of contact rate) respectively. We have 
examined the temporal evolution by displaying successive 
time frame as a movie, but we are unable to represent this 
effect ively on th e printed page. From evolution snapshots 
(Fig. |l(a')]|l(d)[ ), one can see that there is no occurrence 
of spiral waves (fractal fronts) even if the system reaches 
a stable state when the e is in certain interval. The cer- 
tain interval is turned out to be £* < e < £c with the 
parameter value used in the figure and ~ 0.048 and 
Sc ~ 0.305. Here the critical value Sc ~ 0.305 in the spa- 
tial model ([2]), which is more than the value of the local 



dynamic of the chaotic point of system ([T]) [18] (« 0.28). 
This result suggests another possible explanation, which 
is that populations embed and disperse more stabl y in 
space than they do in non-spatial counterparts [23, l25[ . 
As the e {e > Ec-, the domain of chaos) and time in- 
crease, the dynamical patterns with fractal fronts (spirals 
waves) of spatial structures occur, and become larger and 
more stable (see the Fig. [3]). The CA models are generally 
based on qualitative rather than quantitative information 
about the system. It is hardly to detect the density cal- 
culated over the entire lattice when the system is enough 
larger. Hence, in order to investigate quantitatively the 
evolution of the system ([2]), we give the results by one- 
dimensional space. An explicit visualization of spatial 
organization within the lattice is provided by the space- 
time plot of Fig. [2] and Fig. [H 




(a)l step 



(b)5250 steps 




(c)5300 steps 



(d)5350 steps 



FIG. 1: A typical simulation shows four snapshots 
of the evolution in two-dimensional space with the 
parameters /i 0.02, 5 35.84, 7 = 100, /^o 1800 
and £ — 0.23. The figures plot susceptible density 
levels as a space on grey scale. The exposed and 
infected distribution has a qualitatively similar form. 



In Fig. |2(a)| the three time series displaying of the 
density of susceptible, exposed and infected for the first 
10000 steps by CA are given. Self-sustained oscillations 
of the three time series develop (see the Fig. |2(a) ). The 
amplitude of oscillations increase with the increase of 
fluctuating amplitude of the infection rate, e (compared 
the Fig. |2(a)| and 4(a) ). In fact, large oscillations will lead 
to stochastic extinction of the species, when the value of 
fluctuating amplitude is more than in the domain of 
chaos. In Fig. 2(b)|2(c) , spatio-temporal pictures of the 
susceptible and exposed are plotted respectively, where 
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time increases from bottom to top and the horizontal axis 
represents the spatial location. From the Fig |2(b)||2(c) , 
it is clearly seen that the whole system shows the spatial 
period-2 structure when the e is between £* and (later 
we will give the case when e is smaller than a critical 
value £*). 




(a)l step (b)5250 steps 




(c)5300 steps (d)5350 steps 




FIG. 3: A typical simulation shows four snapshots of 
the evolution in two-dimensional space. The parame- 
ters is same as in Fig.dJ /i = 0.02, 5 = 35.84, 7 = 100, 
= 1800 and e = 0.38. The figures plot susceptible 
density levels as a space on grey scale. The exposed 
and infected distribution has a qualitatively similar 
form. 



(a) 



(b) 




(c) 



FIG. 2: The spatial period-2 structure results for the 
system (j2j) in one-dimensional space in cellular au- 
tomaton models and the parameters are the same as 
those in Fig. [H We use a spatial domain of 100 lat- 
tices; 10000 successive time iterations are plotted. In 
(a) we show the species density as a function of time. 
The figure (b) and (c) show the susceptible and ex- 
posed density levels as a function of space and time 
on a grey scale respectively. The behavior of the in- 
fected is qualitatively similar. At this scale the spa- 
tial discretization in not really visible, and therefore 
we have enlarged one region of the plane in figure (c) 
as figure (d), in order to illustrate this discretization. 
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FIG. 4: The spatial period-2 structure results for the 
system (|2]) in one-dimensional in cellular automata 
models and the parameters are the same as those in 
Fig.O We use a spatial domain of 100 lattices; 10000 
successive time iterations are plotted. In (a) we show 
the species density as a function of time. The figure 
(b) and (c) show the susceptible and exposed density 
levels as a function of space and time on a grey scale 
respectively. The behavior of the infected is qualita- 
tively similar. At this scale the spatial discretization 
in not really visible, and therefore we have enlarged 
one region of the plane in figure (c) as figure (d), in 
order to illustrate this discretization. 



To furtherly investigate the impact of the fluctuat- 
ing amplitude on the dynamical patterns with fractal 
fronts (or spiral waves) in two-dimensional and the spa- 
tial period-2 structure in one-dimensional space respec- 
tively, we study the case when the SEIR model is deeply 
in the domain of chaos and out of the domain of chaos. 
The evolution of system ([2|) is shown with e = 0.38 in 
Fig. [3] and Fig. H In Fig. |3(b)||3(d)| three snapshots 
are taken at 5250, 5300 and 5350 steps respectively. In 
Fig. |3(b)| the rotating spirals are not recognizable due to 
the irregular interfaces. However, the spiral formation 
becomes visible when the interfacial roughness grows by 
the infect ed invasion as demonstrated in Fig. |3(d)[ In 
Fig. |3(d)| one can easily identify the vortices and anti- 
vortices rotating clockwise and counterclockwise, respec- 
tively. We have to emphasize that this pattern cannot 
be characterized by a single length unit (e.g., correlation 
length) because the main features of spirals (armlength, 
average curvature, average distance, etc.) depend on the 
model parameters. But it may be analysed by using some 
geometrical features method [26| . The armlength of these 
spiral waves are broad and do not easily break, resulting 
in periodical recurrence of epidemic waves. The spon- 
taneous formation of spiral waves means the regularly 
recurrent infection waves (Fig. [3(iOp(d)] ). Similarly the 



irregular spiral waves can also be observed even when 
the fluctuating amplitude is much more than the critical 
value £c. These results are not shown in this paper. 
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FIG. 5: The spatio-temporal evolution of the system 
([2]) in one-dimensional in cellular automata models 
for the £ < case, £ = 0.035. In (a) we show the 
species density as a function of time. The flgure (b) 
is plot susceptible density as a function of space and 
time on a gray scale. The behavior of the exposed 
and infected is qualitatively similar. 

Figure [5] shows the time evolution of the density of 
the species and the spatio-temporal conflgurations of the 
system (j2j) at £ = 0.035. The spatio-temporal evolution 
of the expose and infected are similar with the suscepti- 
ble's (Figure [5 (b)[ ). It can be clearly noticed that spatial 
period-2 structure disappear and the stationary state is 
a flxed point with the decreasing of £. The situation 
corresponds to a low and persistent endemic infection in 
Fig. 5(a)[ The oscillations decay to the flxed point when 
the £ is smaller than the critical value case. The oscil- 
lations decay because these infection clusters grow, the 
availability of infected hosts per susceptible host is re- 
duced, decreasing the number of new infectors. In this 
case, the spontaneous formation of dynamical patterns is 
qualitatively similar in Fig. [3] in two-dimensional. 



IV. CONCLUSIONS 

A realistic spatial epidemic with the individuals ran- 
domly moving in its neighborhood has been modeled us- 
ing cellular automata. Our simulations demonstrate that 
the recurrent infectious waves exist and persist in an ex- 
haustive spatio-temporal. We have investigated the dy- 
namical patterns of system (|2]) in one and two dimensions 
perspectively. We show that the spiral waves recur peri- 
odically and the recurrence is insensitive to the change of 
the fluctuating amplitude £ within the domain of chaos 
(the fluctuation amplitude £ > £c). Moreover the dynam- 
ical patterns with fractal fronts grow stably. The system 
also shows spatial period-2 structure in one dimension 
when £ is between and £c outside the domain of chaos. 
It is interesting that the spatial period-2 structure will 
break and transform to spatial synchronous conflguration 
in the domain of chaos. Our results conflrm that popula- 
tions embed and disperse more stably in space than they 
do in non-spatial counterparts. 
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